Experimental digital quantum simulation of quantum tunneling 
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It is well known that quantum computers are superior to classical computers in efficiently sim- 
ulating quantum systems. Here we report the first experimental simulation of quantum tunneling 
through potential barriers, a widespread phenomenon of a unique quantum feature, via NMR tech- 
niques. Our experiment is based on the digital simulation algorithm that discretizes continuous 
space and time into lattices, and requires only two spin-1/2 nuclei without the need of ancillary 
qubits. The occurrence of quantum tunneling through a barrier is clearly observed through the 
experimental results. This experiment has clearly demonstrated the viability of quantum simula- 
tions for a wide range of quantum phenomena with not only discrete, but also continuous degrees 
of freedom. 
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Quantum computation has become a subject of intense 
investigation since Feynman discussed the likelihood of 
simulating one quantum system by another [l|. Recent 
years have witnessed fruitful results in the development of 
quantum computation. Not only has Feynman's conjec- 
ture been confirmed, but it has also been demonstrated 
that quantum computers can solve problems, such as 
Shor's quantum algorithm for integer factorization Q 
and Grover's quantum search algorithm p|], with a level 
of efficiency beyond the capability of classical computers. 
As conjectured by Feynman, simulation of the dynamics 
of quantum systems is one of the most important aims 
of quantum computation. Quantum simulations, such 
as the simulation of a many-body interaction Hamilto- 
nian the dynamics of entanglement [7, 8], quantum 

phase transitions M. flO|. and the calculations of molecu- 
lar properties [Tll.ll2j| have been demonstrated in experi- 
ments. Quantum simulations are especially attractive for 
providing an exponential improvement compared to clas- 
sical computations, which makes it a promising candidate 
in the analysis of chemical reactions [13l Il4j. 

As a unique fundamental concept in quantum mechan- 
ics, quantum tunneling plays an essential role in many 
quantum phenomena, such as the tunneling of supercon- 
ducting Cooper pairs [15| . Additionally, its properties 
are widely applied in modern devices and modern ex- 
perimental techniques, such as the tunnel diode [l6[, the 
scanning tunneling microscope [17] and so on. Quan- 
tum tunneling is of continuous interest. Sornborger (l8| 
proposed a small-scale digital simulation algorithm for 
demonstrating a one-particle wave function tunneling in 
a double- well potential using a quantum information pro- 



cessor. In this paper, we realized the algorithm via the 
room temperature liquid state NMR of two qubits. Using 
the digital simulation algorithm, the continuous process 
of one-dimensional tunneling of a particle through a po- 
tential barrier is clearly demonstrated. Our experiment 
has shown that with very few qubits, interesting tunnel- 
ing dynamics is simulated with a gate count that is within 
reach of current quantum architectures. 

In the digital quantum simulation (l8Ti20l|. the one- 
dimensional wave function ip (x, t) of a single particle 
moving in square-like potential with Schroinger equation 
as: 
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where P and X are momentum and position operators, 
respectively. Throughout the text we set H to be 1. In 
the Schroinger picture the evolution of the wave function 
with time can be straightfowardly given as: 



\^(x,t + At)) 



which can be further decomposed according to the Trot- 
ter formula [IlGJl: 

| V 0, * + At)) = \ e - l & At e- lV ^ )At + O (At 2 )] \i/> (x, t)) 

(3) 

For a sufficiently small time step At, Eq. j3| can be 
approximated as: 



IV- (x, t + At)) = e -^ At e~ lVix)At |V (x, *)> 
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Similarly, the continuous degree of coordinate x is also 
discretized. Suppose ip (x, t) is continuous on the re- 
gion < x < L, with a periodic boundary condition 
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if) (x + L,t) = tp (x, t). x is discretized into a lattice with 
spacing Al and the wave function is stored in an n-qubit 
quantum register 



\if>(x,t)) ->■ 1>(xk,t)\k) 



(5) 
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where Xf. — (k + |)Ai, Al — -M;, and |fc) is the lattice 
basis state corresponding to the binary representation of 
the number k. Apparently, Eq. ([S]) gives a good approx- 
imation to the wave function in the limit n — > oo. 

As a small-scale demonstration, a 2-qubit simulation is 
investigated. In this case, the four basis states |00), |01), 
10} and 1 11) register the four lattice points 1, 2, 3 and 4 
as discretized position variables of the particle. 

Next we construct the kinetic and potential operators 
for the discretized wave function. Because the potential 
operator V(X) is a function of the coordinate operator 
X, it is diagonal in the coordinate representation, though 
discretized now. In an n-qubit discretized grid, V, a diag- 
onal matrix in the coordinate representation, can always 
be decomposed as 



V = 



>*fc=l "in 
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where each term in the summation corresponds to a diag- 
onal entry with coefficient c^^...^, 03 is the Pauli matrix 

a z = ( q ^ ), and 174 = I is the identity matrix in two 

dimensions. One simple example of a potential for a 2- 
qubit system with the form of Eq. ([5]) is 



V = V I ® o z 



(7) 



which corresponds to a double-well potential of ampli- 
tude 2Vo, with the two peaks Vb at 00 and 10, and the 
two valleys — Vb at 01 and 11. This double-well potential 
can therefore be implemented using only a single qubit 
gate 
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Similar to Eq. ([5]) , the wave function in the momentum 
representation can be written as 



where 



2"-l 

\4> (?,«))-► E ^fo'Oli) 



s$ j ^ T 



(9) 




is the eigenvalue of the momentem operator P p (the sub- 
script 'p' indicates the momentum representation) corre- 
sponding to the j-th basis The momentum operator 



is expressed in spectral decomposition form as 

p p = E %i \i) 01+ E Sr^-J) li) 01 (10) 
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Therefore the kinetic energy operator in the coordinate 
representation can be written via quantum Fourier trans- 
form (QFT) as 



pi p2 

£_ = F-'-S-i-Z-SF, 
2m 2m 



(11) 



where the QFT is expressed in terms of the bit swap 
gate S, and the bit-swapped quantum Fourier transform 
operater F. 

For the 2-qubit system, Eq. (fTTj) can be specifically 
expressed as 



P 2 tt2 

fL = F ~ 1 s- 1 ^ hsf, 



(12) 
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where we have taken m = 1/2, II = (diag[0, 1, 2, — l]) 2 , 
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and 



F = 5- 1 i ]T e^ fe |j)(fc| 



(13) 
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-F can be readily implemented in quantum circuits via a 
series of Hadamard gates and controlled-phase gates [2lj 



F — H 2 R^Hi, 



(14) 



where Hi and i?2 are Hadamard gates on the first and 
second qubits, respectively, and Rtl = diag [1,1, 
is the two-qubit controlled phase gate. The operator 

■ P 2 A 

e~ l 2^ can be expressed as 



-issrAt = F -i DF _ F^^^ZiZaF, (15) 



where 



r 3^ 2 A , 

<IV = expH— i^At], 



i?,r = RilRjl and 



Zi = exp[i— <T Z ®/Ai], 
3vr 2 

Z 2 = cxp[i— /(gxr^At]. 



(16) 

(17) 
(18) 



Operation exp(— z(27r/4) 2 diag[0, 4, 1, 1] At) is realized 
through ^^ZiZ^, which is also a bit-swapped operator, 
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up to an overall phase factor. Thus the time-dependent 
Schrodhiger Eq. (Q| can be rewritten as 
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i>(x k ,t + At) \k) = F~ l ^Z l Z 2 FQ V V (x k ,t) \k) 
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, (19) 

In the experiment, we studied the time evolution of 
a particle moving in a double-well potential in a digital 
algorithm in a 2-qubit system via the room temperature 
NMR system. A d6-acetone solution of 13 C-labeled chlo- 
roform at 295K is used where 1 H (with a resonance fre- 
quency 400.130 MHz) encodes the first qubit and 13 C 
(with a resonance frequency 100.613 MHz) encodes the 
second qubit. 

From Eq. Q19p the one time step evolution quantum 
circuit for our simulation is straightforwardly obtained, 
as drawn in Fig. [TJ The corresponding pulse sequences 
for implementing the operations in the quantum circuit 
are exhibited in Fig. [5J 
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Figure 1: (Color online) Quantum circuit for one time step 
simulation. |t/>i) and \tjj2) are the input states of thee first and 
the second qubit, respectively. 

In our experiments we set the time interval to be 
At = r-L and the amplitude of the potential 2Vb = 20. 
We simulated the situation in which the particle is ini- 
tially trapped inside one of the two wells by preparing 
the pseudo-pure state 1 01) from the thermal equilibrium 
[22l | as the initial state with an experimentally realized 
high fidelity / = 99.89% (see Fig. [3]) which is given by 



trace(p ex p t h) 



^trace{p 2 ex )trace(p 2 h ) 



(20) 



where pth is the theoretical density matrix, and p ex is the 
experimental density matrix. 

In order to observe the tunneling process of the par- 
ticle, nine experiments have been carried out in which 
the pulse sequence in Fig. [2] has been performed 1 ~ 9 
times. Quantum state tomography (QST) is performed 
on the de nsity matrices of the final states after 1 to 8 
steps ([23L l24l|). The diagonal elements that correspond 
to the simulation of the probability distribution of the 
particle are illustrated in Fig. [5J compared with the the- 
oretical results. The full density matrix of the final state 
after 9 steps in our experiment is shown fully in Fig. [4] 
along with the theoretical results. The fidelity of the 
experimental density is determined by Eq. (|20l) to be 
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Figure 2: (Color online) Pulse sequences to implement the 
circuit in Fig. [TJ where (a), (b), (c), (d) realize the four 
operations F, D, F _1 and Q respectively. The time peri- 
ods il — = 580.9/xs and t2 — -^j = 365.0/is represent 
the free evolution durations under the J coupling (for chlo- 
roform J = 215.2Hz). The bars represent single-qubit rota- 
tion pulses with the different angles as shown to the right, 
and seven different fillings indicate different rotation angles 

as 9± = |, 6>2 = 7T, #3 = f , 84 = -§q, 5 = Jq, 06 = fo> 

B*i = 2, respectively. 
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Figure 3: (Color online). The real and imaginary parts of (a) 
the density matrix of the |01) state; (b) the density matrix 
reconstructed by QST for the initial 1 01} state in our experi- 
ment. 



95.48%. As is depicted in Fig. [6] (b), the whole evolu- 
tion process of a single particle's probability distribution 
in the double-well potential can be obtained through our 
9-step experiment. It can be clearly observed that the 
particle tunnels through the potential barrier between 
the two wells, while its probability to be found in the 
barrier is scarce, which accords well with the theoretical 
calculations in Fig. [5] (a). 

For the sake of comparison we also studied the simula- 
tion of the evolution of a free particle with zero potential 
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(a) Real Imaginary 




Figure 4: (Color online). The real and imaginary parts of (a) 
the theoretically calculated density matrix after nine steps 
of evolution in a double well potential; (b) the experimental 
density matrix reconstructed by QST after the nine steps of 
evolution in a double well potential. 
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Figure 5: (Color online). The real and imaginary parts of (a) 
the theoretically calculated density matrix after nine steps of 
evolution without double well potential; (b) the experimen- 
tal density matrix reconstructed by QST after nine steps of 
evolution without the double well potential. 



using approximately the same experimental schemes and 
parameters, except for the removal of the potential op- 
erator Q in the circuit, and its corresponding RF pulse 
sequences were also performed from one to nine times. 
The final state reconstructed by full QST is shown in 
Fig. [5] with a fidelity of 96.82% according to Eq. ([20]). 
The results of these nine experiments, together with the 
theoretical calculations, are plotted in Fig. [7] to illus- 
trate the evolution of a free particle. It is not surprising 
to find that the probability distributed more evenly on 
all the four lattice points, considering the particle is free. 

Our experiments have simulated the fundamental 
quantum phenomenon of tunneling and reflected an ev- 
idential difference between the two situations with and 
without the double-well potential. If more qubits are 
used in the simulation, more specific dynamic phenom- 
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Figure 6: (Color online). A single particle's probability dis- 
tributions as a function of time for the first nine steps for a 
particle in a double well potential. The two potential wells 
are at |01) (lattice point 2) and [11) (lattice point 4). (a) The 
theoretical predictions; (b) The liquid NMR experimental re- 
sults. The four bars (from left to right) at each step indicates 
the particle's probability on the state 1 00) (blue), 1 01) (green), 
1 10) (yellow), 1 1 1) (red), from left to right. The initial state 
is prepared in |01). As is illustrated in (a) and (b), the par- 
ticle tunnels from |01) to |11). Our experiments agree quite 
well with theoretical results depicted, clearly manifesting the 
tunneling phenomenon. 



ena of quantum tunneling can be observed, for example, 
with higher position resolution in wells or in barriers, 
the oscillation of the state in wells or in barriers can be 
observed. An TV-lattice-point simulation can be realized 
with only log2N qubits, which means the scheme is effi- 
cient and easily implemented in experiments. The exper- 
imental realization of multi-qubit simulation of quantum 
tunneling will be studied in our future work. 

It should be emphasized here that although the real 
evolution of the particle takes place in a continuous space 
with infinite dimensions, the quantum computer, via only 
a few qubits with limited dimensions (in our case four di- 
mensions), is already capable of undertaking some basic 
yet fundamental simulation tasks, such as quantum tun- 
neling. The result has revealed the amazing power hidden 
under the qubits and the promising future of quantum 
computations. 

In summary, we accomplished a small-scale demonstra- 
tion of the quantum tunneling process on the 2-qubit 
NMR system based on the digital quantum simulation 
method. To our knowledge this is the first demonstra- 
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tion work on quantum tunneling simulations via an NMR 
quantum information processor. The experimental re- 
sults and the theoretical predictions are in agreement. 
The overlaps of the final states after nine steps of evo- 
lution (with and without the double-well potential) be- 
tween the experimental results and the theoretical pre- 
dictions are about 95%, which manifests the high quality 
of the experimental results. 
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Figure 7: (Color online). The free particle's probability dis- 
tributions as a function of time for the first nine steps, (a) 
The theoretical predictions; (b) The liquid NMR experimen- 
tal results. The four bars at each step indicates the particle's 
probability on the state 1 00) (blue), 1 01) (green), 1 10) (yellow), 
1 1 1 ) (red), from left to right. The initial state is prepared to 
be 1 01) . As is illustrated in (a) and (b), the particle's prob- 
ability becomes more evenly distributed among all the four 
states. Our experiments agree quite well with theoretical re- 
sults depicted. 
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